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1. Introduction 

This is a talk based on [1]. The action of the 2-d 0(3) model in the continuum is defined as 



where g 2 is the coupling and the unit-vector e(x) £ S 2 is the field variable. The 2-d 0(3) model 
shares many features with QCD: Both have important properties including asymptotic freedom and 
a non-perturbatively generated massgap. One can define a topological charge [2] 



which is an integer number and describes how many times the sphere S 2 gets covered by e(x), when 
x runs over a periodic space-time. There are instantons, which are field configurations with a non- 
zero topological charge Q[e\, which minimize the action S[e\. Because of the presence of non-trivial 
topological sectors we can add the 0-vacuum term i6Q[e\ to the action. Here, we show that the 
parameter 6 is indeed a relevant parameter in the action, which means it does not get renormalized 
non-perturbatively. We will show that at each value of 6 there is a different continuum theory. 

Haldane showed in the large spin S limit that the low-energy effective field theory of a (l+l)-d 
spin chain is the 2-d 0(3) model. He further conjectured that integer spin chains correspond to a 
vacuum angle 6 = 0, while half-integer spin chains correspond to 6 = n [3]. This conjecture has 
been confirmed numerically for 6 = in [4] and within statistical errors also for 6 = 71 [5]. For 
spin chains it is well known that for half integer spins the mass gap vanishes [6]. 

The 2-d 0(3) model has been studied analytically by the exact S-matrix theory. For 6=0, this 
theory predicts a mass gap m{6 = 0) = -A^, which is generated non-perturbatively and has been 
verified by lattice simulations. Furthermore it is conjectured that the low-energy effective field 
theory of the 2-d 0(3) model at 6 = K is the WZNW model [7]. The exact S-matrix theory even 
solves analytically the 2-d 0(3) model also in the finite volume and provides a prediction for the 
step scaling function [8]. Therefore, it is not necessary to take the thermodynamic limit to compare 
analytic and numerical results. The step scaling function a(s, uq) describes the scaling behavior of 
the renormalized coupling «o = m{L)L in the finite volume with spacial extend L. Here, we will use 
different lattice actions to confirm the exact S-matrix results. Besides the standard action, we use a 
topological action [9], which constrains the maximal angle between neighboring spins. To reduce 
the cut-off effects, we will combine these two actions in the new "constraint action", which has 
also been studied intensively in [10]. In this paper, we also use a variant of a method developed by 
Hasenbusch [11, 12] to simulate 0-vacuum effects in the 2-d 0(3) model with unprecedented per 
mill level precision. For the first time, this numerically confirms the conjectured exact S-matrix of 
the 2-d 0(3) model at 6 = 71 [13] beyond any reasonable doubt, which also implies that the model 
indeed reduces to the WZNW model at low energies. 

On the lattice, field configurations with non-zero topological charge that minimize the action 
are known as dislocations. When the dislocation action is less than a critical value, a semi-classical 
argument suggests that the topological susceptibility %t = (Q 2 ) /V should suffer from an ultra- 
violet power-law divergence in the continuum limit [14]. As we will see later, for each value of 6 
the dislocation problem does not prevent 6 to be physically relevant, which suggests that % t only 
diverges logarithmically. 




(1.1) 




(1.2) 
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2. Lattice Setup 

2.1 The Lattice Actions 

The simplest lattice discretization of the action (1.1) leads us to 

Standard [e\ = 3 Hi ^ ' ^ ' ^ 1 ) 

8 (xy) 

where e x G 5 is the field variable, which is now denned on every site x on the lattice. Besides this 
standard action we also consider the topological action introduced in [9] 



r-n r-i v ^ JO for e x -e y >cos<5 

^topological^ = y,s{e x ,e y ), s{e x ,e y ) = < . (2.2) 

1 00 else 



Here, the maximally allowed angle 5 plays the role of the coupling. The topological action only 
allows configurations where the angle between the field variables on neighboring sites x and y is 
smaller than 8 (e x • e y > cos 8). Otherwise the action is infinite, which means that the corresponding 
configurations are not allowed. All allowed configurations have the same action ^topological [e\ = 
0. As a consequence, this lattice model does not have the correct classical continuum limit, it 
violates the Schwarz inequality between action and topological charge, and it cannot be treated 
in perturbation theory. Despite these various deficiencies this action still has the correct quantum 
continuum limit [9]. 

As it will be shown later, the standard action approaches the continuum limit of the step scaling 
function from above while the topological action approaches it from below. Therefore we combine 
these two actions in order to reduce the cut-off effects. This we do with the an optimized constraint 
action 



l ^ 

o2 



e x ■ e y for e x ■ e y > cos 8 



^constraint!?] = Y. s'(e x ,e v ), s' (e x ,e y ) = < g2 x y * y _ (2.3) 

S I 00 else 

Here, g 2 is again the coupling, while 8 is a fixed parameter that is tuned to a value which minimizes 
the cut-off effects. For 8 = % this optimized constraint action reduces to the standard action (2.1). 
On the other hand, if we send g 2 — > oo, we obtain the topological action (2.2). 

To implement the topological charge Q[e\ in a discretized form, we triangulate the lattice as 
shown in figure 1. Then each triangle (xyz) gets mapped to an oriented area A /^, z \ = 47iq^ xy ^ on the 
sphere S 2 defined by e x , e y , and e z . If we sum up the areas qi^d of ai l triangles, taking in account 
their orientation, we obtain the topological charge 

Q[A= E?W- (2-4) 

(xyz) 

Using this definition, the topological charge is an integer number, as long as we sum over all 
triangles in a periodic lattice. The area of one single triangle (xyz) can be calculated as 

Re 2ni q{xyz) = l+?x . ?y + Sy .e z + e z - e x + ie x ■ (e y xe z ), R>0. (2.5) 



As before, we can add the 0-vacuum term i6Q[e\ to the action, where Q[e\ is now the dis- 
cretized topological charge. 
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A(t u t 2 ) 



Figure 1: Triangulated square lattice: the triangles (xyz) in the shaded area A(t\ ^2) carry the topological 
term idq^ xyz y 



2.2 Observables 

To compare results with the exact S-matrix theory, we calculate the step scaling function 
a(s, t<o) [8], where s is a rescaling factor. Starting on a volume with spacial extent L, one ad- 
justs the coupling g 2 in order to obtain the renormalized coupling uq = m(L)L, where m(L) is the 
inverse of the finite volume correlation length m(L) = l/<§ (L). Then, one measures the renormal- 
ized coupling ui = m(L')L' on the scaled volume with spacial extend L' = s-L. Setting the rescaling 
factor s = 2 means that we measure the quantity 



a(2, u) =2m(2L)L 



m(L)L=UQ 



(2.6) 



on a volume which has twice the spacial extent L' = 2L than the original. To measure the finite 
volume correlation length (L), we calculate the correlation function 

C(h,t2-,0) = z{h ) 2 . d) (jlj^) E( tl ).E(t 2 )xe- s ^ w ^ ~ e-" 1 ^-"), (2.7) 

where Q(ti,t 2 ) is the (in general non-integer) topological charge between time-slice t\ and t 2 , 
m(Q,L) is the inverse of the finite volume correlation length m(6,L) = 1/<§(0,L), E(t) is the 
average spin in timeslice t and Z(t\,t 2 ; Q) is the partition function 

^( ? ) = ^rL<W Z(t l ,t 2 \e)= (Uj s2 ^ exp(-S[e]+idQ{h,t 2 )). (2.8) 

To sample Monte Carlo configurations with a positive weight, we need to absorb the phase factor 
in the observable. This reweighting is obtained as follows 



x E( h ).E(t 2 )) «p(»e6fti,fa)); fl=0 c(fl>f2> o )z(fl>f2> fl) /z(0 ) 

u ' 2 ' j " <ex P (/ee(^ 2 ))> e=0 z( tl ,t 2 ,e)/z(o) ■ 

The cut-off behavior of the lattice step scaling function l,(8,2,uo,a/L) is known for 6 =0[15]: 

a 2 

£(0=0, 2, u ,a/L) = o(0 = 0, 2, u ) + [filog 3 (Z» + Clog 2 (Z» + ...], (2.10) 
where B and C are fitting parameters and L/ais the number of lattice points in the spatial direction. 
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To obtain high accuracy data, we use an improved estimator first proposed by Hasenbusch [11], 
which was further developed in [12]. Instead of measuring the estimator (&(ti) -E(t 2 )\ e ie Qfa< t2 \ 
we integrate out a certain rotation symmetry. This can be done by rotating all spins after timeslice 
t with a rotation matrix X(t) : euf) 1 — > X(t)e(x.t') f° r a ^ ? ' > f - Using open boundary conditions 
in the temporal direction, the action S[e\ and the topological charge Q[e\ only change in this time 
slice t. For each timeslice t\ < t < t 2 we now choose four different rotation matrices Xo(t) = 
1, X\(t), X 2 (t), Xi(t), which all do not change the action S[e] (for the construction of see 
[12]). So in each timeslice t\ < t < t 2 , we insert a matrix 

X{t)=YXi{t) ex V (ieQ(t,t + l,Xi(t))), (2.11) 

i 

where Q(t,t + is the topological charge after all the spins et x j') for t' >t have been rotated 

with the rotation matrix With this, we can write improved estimators as 

(E(h)-E(t 2 ^ exp(i6Q(t h t 2 )) ~ E{h) X{h + 1) X{h +2) • • -X(t 2 ) E{t 2 ) (2.12) 

esp(iOQ(t u t 2 )) ~ ... (j^ e ieQ(t2-ht 2 Mt2))\ 

i i 

3. Results 

By measuring the step scaling function on different lattice sizes L/a, we can plot the cut-off 
behavior of the different lattice actions as it is done in figure 2. The curves for the standard action 
and for the topological action are fits of the function (2. 10), where a, A and B are fitting parameters. 

For 8=0 the continuum limit (horizontal line) is predicted by the exact S -matrix theory 
to be at a(0,2, 1.0595) = 1.2612103 [16]. As one can see, the standard action approaches this 
continuum limit from above, while the topological action approaches the continuum limit from 
below. The maximum angle 8 of the constraint action is tuned to the continuum value at a lattice 
size L/a = 10 and takes an optimal value at cos(5) = —0.345. Even for finer lattices the cut-off 
effects of the constraint action are remarkably small and only at the per mill level. 

At 6 = n the exact S-matrix theory again predicts a continuum result (horizontal line) of 
o(7l,2, 1.0595) = 1.231064 [17]. This result is confirmed by extrapolating the standard action or 
the constraint action. For the constraint action, we still use the maximal angle 8, which has been 
tuned at 6 = 0. Also in this case the cut-off effects of the constraint action are extremely small. 

For 6 = 71 /2, we do not have an analytic prediction of the continuum limit and we do not know 
whether the function (2.10) is still applicable to describe the cut-off effects. Nevertheless, we fit this 
function to the cut-off effects of the standard action, which gives a small # 2 /d.o.f. Extrapolating 
to the continuum limit agrees with an estimator of the continuum limit for the constraint action 
(horizontal line). Again, the constraint action has a maximum angle constraint 8, which has been 
tuned at 6 = 0, but still shows an extremely good cut-off behavior. 

Figure 2 also shows that all continuum limits (at 6 = 0,7t/2,7t) are significantly different, 
which means that each value of 6 indeed corresponds to a different theory. Hence 6 is a relevant 
parameter that does not renormalize to or 71 non-perturbatively, as one might have expected due 
to the presence of dislocations. This can also be seen in figure 3, where we plot m(8,L)L as a 
function of 6, keeping the coupling fixed. For 8 = 0, we fix the value to m(0,L)L = 1.0595, at 
d = % one obtains m(%,L)L = 1.048175 [17]. 
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Figure 2: Cut-off dependence of the step scaling function 1.(9, 2, uq = 1.0595, a/L) for three different lattice 
actions: the standard action, the topological lattice action of [9], and the optimized constraint action with 
cos 8 = —0.345, as well as for three different values of 9 = (top) and 9 = f , n (bottom). The lines are fits 
based on eq. The horizontal lines represent the analytic continuum results for 9 = [16] and 9 = n [17], 
and the fitted continuum value for 9 = 5 . 



4. Conclusion 

Using an efficient modification of Hasenbusch's improved estimator allowed us to show that 
the 6 parameter in the 2-d 0(3) model is indeed a relevant parameter and does not get renormalized 
non-perturbatively. We found a different continuum limit for each value of 6 (here shown for 
6 = 0,7l/2,7r). Dislocations do not spoil this continuum limit even for 6^0. We confirmed the 
exact S-matrix conjecture for the step scaling function at 6 = %, which also implies that the model 
indeed reduces to the WZNW model at low energies. 

In the simulations we used different lattice actions. The constraint action combines the stan- 
dard action with a topological action. The topological action restricts the maximal angle between 
neighboring spins and is thus invariant under small deformations of the field. In the constraint 
action this maximal angle is tuned in order to reduce the cut-off effects. Our results show that only 
a single tuning is necessary (at L/a = 10, 6 = 0) to once and for all fix this maximal angle. Other 
lattice sizes or other values of 6 have extremely small cut-off effects of the step scaling function 
which all lie in the per mill level or beyond. 
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Figure 3: The ^-dependent massgap Lm(9,L) at Lm(Q,L) = 1.0595 using the optimized constraint action 
for L = 24a, compared to the analytic result at 9 = n [17] (cross). 

We are indebted to J. Balog for providing exact values of the step scaling function at 6 = 
71 prior to publication. This work has been supported by the Regione Lombardia and CILEA 
Consortium through a LISA 201 1 grant, as well as by the Schweizerischer Nationalfonds (SNF). 
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